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O [ Abstract 

■ The net-baryon density at midrapidity is proposed as an order 

^ , parameter in the search for the QCD critical point in heavy ion col- 

| lisions. As a function of the initial energy and the total number of 

participants, this quantity obeys a scaling law, dictated by the critical 
exponents of the appropriate universality class. The corresponding 
scaling variable specifies the proximity of a given experiment to the 
jJJ ■ critical point. Within this framework, measurements at the SPS are 
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discussed and predictions for RHIC and LHC are extracted. 
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1 Introduction 



A remarkable property of the QCD phase diagram (p — T) is the existence 
of an endpoint along the critical line of the first order quark-hadron phase 
transition [|TJ . It defines a critical point of second order, belonging to the uni- 
versality class of a 3d Ising system and located on a line of nonzero baryonic 
density p = p c . The QCD critical point is associated with the chiral phase 
transition in the sense that it is the remnant of a tricritical point correspond- 
ing to the chiral limit m u = uid = 00. In other words, the existence of 
a second-order critical point, at nonzero baryonic density, is a fundamental 
property of real QCD with small but nonzero quark masses (m u ,m^ 7^ 0). 
The QCD critical point communicates with the hadronic world through the 
fluctuations of a scalar field (cr-field) which carries the quantum numbers of 
an isoscalar (a- meson) as the manifestation of a quark condensate, a ~ (qq) , 
in thermal environment. At the critical temperature (and infinite volume) 
the isoscalar has zero mass, in order to provide the infinite wavelength mode 
required by the divergence of the correlation length. It remains, therefore, 
stable in the thermal environment of a nuclear collision as long as the tem- 
perature stays close to the critical value. At the freeze-out stage the mass 
of the ex-field may reach the two-pion threshold (m a ~ 2m n ) and become 
accessible to observation |2|, |3|. In the effective theory of the QCD critical 
point, the classical cr-field is a natural order parameter, the fluctuations of 
which obey scaling laws dictated by the critical exponents of the 3d Ising 
system (77 « 0, (3 « ~, 8 « 5, v ~ |). In a baryonic environment, how- 
ever, the chiral condensate is expected to have, at T = T c , a strong depen- 
dence on the net-baryon density, driving the cr-field close to zero for p « p c : 
{q.q)p ~ ^ (99)o + 0[(P ~ Pc) 2 ] where A is a dimensionless constant 

of the order of unity [|J. This dependence suggests a new order parame- 
ter, m = p — p c , associated with the critical properties of the baryonic fluid 
created in a quark-hadron phase transition. In fact, approaching the critical 
point in the phase diagram, both the cr-field fluctuations and the fluctuations 
of the order parameter m(x) obey the same scaling laws ((qq) p ~ m(x)). 

In this work we exploit the scaling properties of the order parameter 
m(x), properly adjusted to measurable quantities, in heavy ion collisions. 
For this purpose we consider net-baryon production in collisions of heavy 
nuclei with total number of participants A t and initial energy corresponding 
to a total size in rapidity Ay = L. The paper is organized as follows: In 
Section 2 we discuss in detail the scaling properties of the net-baryon fluid 
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near the critical point. In Section 3 we present shortly phenomenological 
consequences of the presence of critical fluctuations in the baryonic sector. 
Finally in Section 4 we present our conclusions concerning the possibility to 
approach and observe the QCD critical point in heavy-ion colliders. 

2 Net-baryon scaling 

The created baryons in the process of quark-hadron phase transition occupy 

.1/3 

a cylindrical volume with transverse radius R± ~ AJ_ and longitudinal size 
L (in rapidity). The parameter A± specifies the effective number of par- 
ticipants, contributing to the transverse geometry of the collision, and it 
is assumed A± « valid both for central (Aj_ = A m i n ) and non-central 
collisions. Projecting out the net-baryon system onto the longitudinal di- 
rection we end up with a Id liquid confined in a finite rapidity region of 
size L with local density p(y) = J^i^ , directly related to the measurable 

net-baryon density in rapidity rib(y) = Putting R± = RoA 1 / 3 , we in- 
troduce a characteristic volume V Q = TtR 2 Tf in terms of the freeze-out time 
scale Tf > 6 — 8 fm || and the nuclear-size scale R Q which contains also 
any growth effects near the critical point (R > 1.2 fm). ). Using V' 1 as a 
scale for baryonic, freeze-out densities, the order parameter m(y) of the Id 
baryonic liquid is written: 

m(y) = Al 2/3 n b (y) - p c , 0<y<L (1) 

In what follows, our basic assumption is that the deconfined phase of quark 
matter, in heavy-ion collisions, approaches the QCD critical point in local 
thermal equilibrium. In the framework of inside-outside mechanism, this 
process is implemented by considering the isothermal space-time hyperbo- 
las t 2 — x| = t 2 and employing the corresponding rapidity variable as the 
appropriate longitudinal coordinate. In this description, the assumed local 
equlibrium in the conventional geometry (xtt,x±,t) is translated as global 
equilibrium in the new geometry (y, x±, r) and one may consistently impose, 
near criticality, static scaling laws on the order parameter m(y). In particu- 
lar, the expected singularity aXT = T C is incorporated in a general expansion 
of the form: 

m(y) = i? [F (y/L) + tF 1 (y/L) + ...] (2) 
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where t = and (3 is the appropriate critical exponent (/3 ~ |). The 

leading term F Q (y/L) in the expansion (Q) has a universal power-law be- 
haviour near the walls (y = 0, L), shared by all fluids belonging to the same 
universality class and confined in a finite, one-dimensional region |], (7|. At 
midrapidity (y ~ |f), eq.(H) gives the deviation of the measurable, bulk den- 
sity of net baryons, from the critical value p c , as we approach the critical 
point, not along the critical isochore (p = p c ) but along the freeze-out line: 

A- ± 2/3 n b = Pc + tf[F + t f F 1 + ...} (3) 

where tf — C T i , n b = n b (L/2), Fi = Fj(l/2) and Tf is the freeze-out 
temperature. Integrating now eq.(^) m the interval < y < L we obtain at 
T = T f 

Al 2/z A t L- 1 = Pc + t%I +t f h + ...) (4) 

where = fo Fi(£)d£. Introducing the variable z c = A^ 2 ^ A t L~ l we find 
from eqs.(^) and (§) a scaling law for the net-baryon density at midrapidity: 

A^_ 2/3 n b = ^(z c ,p c ) ; z c >p c (5) 

where the scaling function *$?(z c , p c ) has the property *$?(z c = p c , p c ) = p c . In 
the crossover regime z c < p c , where critical fluctuations disappear, the local 
density at midrapidity is, to a good approximation, n b w A t L~ l suggesting a 
continuous extension of the scaling law (|^) in this region with ^(z c , p c ) = z c 
(z c < p c ). It is of interest to note that although the scaling function ^f(z c , p c ) 
is continuous at the critical point (z c = p c ), the first derivative is expected 
to be discontinuous at this point in accordance with the nature of the phase 
transition (critical point of second order). 

In order to study in detail the critical behaviour of the baryonic fluid 
in terms of the new variables n b , z c a further investigation of the structure 
of the scaling function *$?(z c , p c ) for z c > p c is necessary. For this purpose 
we approximate the nearby part of the freeze-out line, close to the critical 
point (tf <^ 1), by truncating the series @ keeping only the next to the 
leading term F^y/L). As a result, the equations (^) and (f|) are written 
correspondingly: 

Al 2/3 n b w pc + tfiFc + ifFj (6-a) 

Zc-Pc ~ Mlo + tfh) (6-b) 
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Neglecting terms of order 0(tV) in the power expansion of the quantity (I + 
tfliY'P and combining eqs.(6) we finally obtain: 

A- ± 2/3 n b = Pc + ^[f( Zc ,p c )f + C[f(z c ,p c )f +1 (7-a) 
■fx 



CT 

Al 2/ \ b = Pe + F ^{l + t f —^—) (7-b) 

-f o 



/(2c, Pc) = ^(-! + y/l + 2G(z c - Pc y/P) ; z c > p c (7-c) 



where C = G = -^h-. 

The component F (y/L) which dominates the order parameter m(y) in 
the limit T — > T c is approximately constant in the central region « 
■| , L ^> 1), far from the walls (at the points y = 0,L) due to the ap- 
proximate translational invariance of the finite system in this region. On the 
other hand, approaching the walls, F (y/L) describes the density correlation 
with the endpoints and obeys a universal power law. In summary: 

F o (0 ~ const. (£«-) (8-a) 

F o (0 ~ r P/ " (C>0) (8-b) 
F o (0 ~ {l-t)-W» (£<1) (8-c) 

The solution which fulfils these requirements is: 

fo(0 = <?[£(!- O]^ (0 < e < 1) (9-a) 

F API" „, . T(m)T(n) , . 

— — = H a- ; E n, m) = ^ ; V , 9-b 

Inserting the constant j 2 - into the eq.( [T-b| ), the final form of the scaling 
function ty(z c ,p c ) is obtained: 

^(z c , Pc ) = p c + —j -— g- [/(z c ,p c )] 

+C7[/(z c ,p c )] /J+1 ; ^ c >p c (10-a) 
*(*«Pc) = z c ! < Pc (10-b) 



where f(z c , p c ) is given by eq.(|7-c). It is straightforward to show that the dis- 
continuity of the first derivative of ^(z c , p c ) at z c = p c is a nonzero universal 
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constant: disc Qff-) ~ 1 — f (f ~ §) as expected from the characteris- 
tic properties of a second-order phase transition. Combining eqs.(|^-10) one 
may propose a framework for the treatment of certain phenomenological as- 
pects of the QCD critical point. The scaling law ([|) combined with eq.(10) 
involves three nonuniversal parameters: (a) the critical density p c and (b) 
the constants C, G which give a measure of the nonleading effects, allowing 
to accomodate in the scaling function (10) processes not very close to the 
critical point. We have used measurements at the SPS in order to fix these 
parameters on the basis of eqs.(^) and (10). More specifically, in a series 
of experiments (Pb+Pb, S+Au, S+Ag, S+S) with central and noncentral 
(Pb+Pb) collisions at the SPS || ||] net baryons have been measured at 
midrapidity whereas the scaling variable z c = A~^ A t L~ x y associated with 
these experiments, covers a sufficiently wide range of values (1 < z c < 2) 
allowing for a best fit solution. The outcome of the fit is consistent with the 
choice G ~ and the equations (7) are simplified as follows: 



a —2/3 

A, n b 



Al 2, \ b 

1-2/3 



Pc + -(z c - p c ) + c(z c - p c y 

7T 



p c ^ 1 (i 

7T 1 



A ± n b (y) = p c + 



(z c - p c )L 



71 



+ 0[(z c - p c 



(11-a) 

;n-b) 

(ll-c) 



In eqs.(ll) we have used the approximate values of the critical exponents 
/3 ~ |, ^ ~ | in the 3c? Ising universality class flT0| . We have also added 
eq.( |ll-c| ) which gives the universal behaviour of the net-baryon density rib(y), 
in the vicinity of the critical point (z c — p c 1). The fitted values of the 
parameters in eq.( |ll-b ) are p c = 0.81, C = 0.68 and the overall behaviour of 



the solution is shown in Fig. 1. In turns out that the critical density is rather 
small, compared to the normal nuclear density p Q ~ 0.17 frrr z : p c < 4f 
(R Q > 1.2 fm, Tf > 6 fm), suggesting that the critical temperature remains 
close to the value T c w 140 MeV obtained in studies of QCD on the lattice 
at zero chemical potential 



111. The difference d r 



\z c — Pc| in Fig. 1 is a 
measure of the proximity of a given experiment (A t , L) to the critical point. 
We observe that the central Pb+Pb collisions at the SPS (d c ~ 1.1) drive 
the system into the most distant freeze-out area, from the critical point, as 
compared to other processes at the same energy. In fact, the most suitable 
experiments to bring quark matter close to the critical point at the SPS are: 



5 



S+S (4 « 0.18), S,+Si (4 « 0.20) and C+C (d c w 0.06) ), central collisions. 
One may even reach the critical point (d c ~ 0) at the SPS with medium- 
size nuclei, either in noncentral collisions: S+S, Sj+S; (A t ~ 29) or at lower 
energies: C+C (P pa 130 A^-). Also, the experiments at RHIC, with 
central Au+Au (d c ~ 0.27) collisions, come close to the critical point and may 
even reach it exactly (d c pa 0) either with lighter nuclei or with noncentral 
Au+Au collisions (A t pa 165). Finally at the LHC (^s ss 5.5 Tel/), Pb+Pb 
collisions are expected to drive quark matter into the crossover area (z c < p c , 
d c pa 0.13) where no critical fluctuations occur. At lower energies however 
[y/s pa 1.4 TeV), Pb+Pb central collisions may reach the critical point even 
at the LHC. Obviously, in order to have a clear answer on the proximity 
of a given experiment to the critical point, precision measurements of the 
net-baryon density at midrapidity, in a wide range of energies (between SPS 
and RHIC), are needed, both for central and noncentral collisions. With 
such measurements, the scaling law (§) may be tested in its full extent and a 
sharp determination of the critical density can be achieved through eq.(7a). 

The freeze-out line (n^ versus Tf) given by eq.(llb) is shown in Fig. 2. 
We have used the freeze-out temperatures for the processes f[2| : Pb+Pb 



(T f pa 115 MeV) and S+S (2) pa 145 MeV) at the SPS, in order to de- 
termine the critical temperature, T c pa 145 MeV, and fix the parameter 
I pa 1.86. Again, we notice that the relatively small value of the freeze-out 
temperature in Pb+Pb collisions (at the SPS) is associated with the fact 
that in this process the system of quark matter is not driven close to the 
critical point. On the contrary the freeze-out temperature in S+S collisions 
remains practically equal to the critical value, due to the proximity of this 
process to the critical point, at the freeze-out stage (Fig. 1). Finally, in 
Fig. 3, the rapidity distribution of net baryons in Au+Au collisions at RHIC 
(L pa 11, z c pa 1.1) is shown, relying upon the fact that in this process the 
distance from the critical point is small (d c pa 0.27) and eq. (|ll-cp is valid. 



3 Baryonic critical fluctuations 

Once the phase of quark matter has reached the critical point in a particular 
class of experiments, as discussed in the previous section, strong critical 
fluctuations are expected to form intermittency patterns both in the pion 
and net-baryon sector. As already mentioned in the introduction, the origin 
of these fluctuations can be traced in the presence, at T = T c , of a zero mass 
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field with a classical profile (a-field) which, under the assumption of a phase 
transition in local thermal equilibrium, is described by an effective action in 
3 — d, the projection of which onto rapidity space is written as follows |3|: 



ttR 2 



Ca JSy 



dy 



Kl) + 2c ^ v+l 



Ca = Pc = T- 1 (12) 

Pc 



Equation (12) gives the free energy of the cr-field within a cluster of size 5y in 
rapidity and R± in transverse space. The critical fluctuations generated by 
(|T2| ) in the pion sector have been studied extensively in our previous work , 
therefore, in what follows, we are going to discuss the fluctuations induced 
by the a-field in the net-baryon sector, noting that a direct measurement 
of these fluctuations may become feasible in current and future heavy-ion 
experiments. For this purpose we introduce in eq.(|T2"D the order parameter 
m(y) through the following equations: 



a{y) « Ff3 2 c m{y) ; F 

r -, / o - \ 2 

r c ~ gi / dy 



1 / dm 

2(i»J +92|m 



Pc 
5+1 



; (qq)o ~ -3 fm 3 

; m(y) = f3 3 c m(y) (13) 



where: g x = F 2 (g^r) , 92 = 2C^F 4 . The partition function Z = J V[m]e- r ^ 
for each cluster is saturated by instanton-like configurations [13| which for 
$y — $c l ea d to self-similar structures, characterized by a pair-correlation 
function of the form: 



/ w w , » 5T 1/3 UR\C A 
(m(y)m(O)) « — ^-^ 



F l y s+i 



(14) 



2/3 



The size, in rapidity, of these fractal clusters is 5 C ~ ( ig/pc' 2 J according 

to the geometrical description of the critical systems Integrating eq.(|T4D 
we find the fluctuation (Sn^) of the net-baryon multiplicity with respect to 
the critical occupation number within each cluster, as follows: 



(Sn b ) 



RICA 2 2 / 3 r(i/3) 5/6 

2(3! ) r(l/6) c 



(15) 
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The dimensionless parameter F is of the order 10 2 and the size S c , on general 
grounds (R± ~ 2r c ) is of the order of one (S c ~ 1). This is in agreement with 
the fact that the rapidity separation of two causally correlated space-time 
events is logarithmicaly bounded, 5y < In ( 1 + ^rM , leading, practically, to 

a direct correlation length of the order of one (5y ~ 1). However, in reality, 
the global baryonic system develops fluctuations at all scales in rapidity since 
the direct correlation (0(5 C )) propagates along the entire system through the 
cooperation of many self-similar clusters of relatively small size (5 C ~ 0.35 
and (Srib) m 140). We have quantified this mechanism in a Monte-Carlo sim- 
ulation for the conditions of the experiments at RHIC (in RHIC the size of 
the system is L « 11) in order to generate baryons with critical fluctuations. 
The distribution of the order parameter |m(y)| describing the fluctuations 
of the "critical" baryons in the rapidity space for a typical event as well as 
the corresponding intermittency analysis in terms of factorial moments are 
presented in Fig. 4. The intermittency exponent of the second moment F 2 
in rapidity is found to be s 2 ~ 0.18 which is very close to the theoretically 
expected value (|) of a monofractal 1 — d set with fractal dimension |. The 
quantitative details of this new class of phenomena, concerning critical fluc- 
tuations in the net-baryon sector, and in particular the intermittency pattern 
of the observable net-baryon density, go beyond the purpose of this Letter 
and will be presented elsewhere fi~4| . 



4 Concluding remarks 

In conclusion, we have shown that the baryonic sector in an experiment with 
heavy ions, possess valuable information regarding the proximity and obser- 
vation of the second order QCD critical point. Complementary suggestions 
for the significance of net-baryon fluctuations in order to trace the critical 
line of first order in the phase diagram are described in [[UJ. In particular 
the measurements of net-baryon spectra in rapidity provide a valuable set of 
obsevables in heavy ion experiments, in connection with the phenomenology 
of the QCD critical point. The trend of the existing data at the SPS suggests 
the presence of a critical point of second order in the phase diagram (Figs. 1,2) 
specified by the critical values of temperature and density: T c rs 145 MeV, 
Pc < A scaling law for the net-baryon density at midrapidity nt,, as a 
function of the initial energy and the number of participants, has been es- 
tablished in the neighbourhood of the critical point. The scaling function 
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incorporates the indices of the universality class (critical exponents) whereas 
the scaling variable gives a measure of the proximity of a given experiment 
to the critical point. On the basis of this investigation, the experiments at 
RHIC are very likely to reach the QCD critical point and as a first sign of 
this new phenomenon we have predicted a rather unconventional profile of 
net-baryon density in rapidity (in Au+Au collisions), associated with the 
presence of a critical point, nearby (Fig. 3). Finally we have indicated the 
presence of strong intermittency effects in the rapidity spectrum of the net 
baryons which can in principle be used to reveal experimentally the QCD 
critical point and confirm its universality class (Fig. 4). 
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Figure Captions 



FIG. 1. The scaling law (10) is illustrated together with measurements at the 
SPS. The critical point and the corresponding break in the slope of ^(z c ,p c ) 
are also shown. 

FIG. 2. The freeze-out line, function oft/, is illustrated together with 

measurements at the SPS (S+S, Pb+Pb). The extracted freeze-out points 
for S+Au and S+Ag are also shown in this diagram. 

FIG. 3. The net-baryon profile in rapidity for Au+Au collisions at RHIC 
energies is shown, as predicted by eq. flll-cp . 

FIG. 4. (a) The distribution of |77i(y)| in rapidity for a MC-generated event, 
(b) The first three factorial moments for the event shown in (a) in a log-log 
plot. A linear fit determining the slope S2 (~ 0.18) of the second moment is 
also shown. 
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